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The phase transitions that occur in an infinite-range-interaction Ising ferromagnet in 
the presence of a double- Gaussian random magnetic field are analyzed. Such random 
fields are defined as a superposition of two Gaussian distributions, presenting the 
same width a. It is argued that this distribution is more appropriate for a theoretical 
description of real systems than its simpler particular cases, i.e., the bimodal {a = 0) 
and single Gaussian distributions. It is shown that a low-temperature first-order 
phase transition may be destructed for increasing values of cr, similarly to what 
happens in the compound FexMgi_xCl2, whose finite-temperature first-order phase 
transition is presumably destructed by an increase in the field randomness. 
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I. INTRODUCTION 

The theory of continuous phase transitions is very well established nowadays 
these transitions are characterized by a divergence of the correlation length, in such a way 
that some microscopic details of the system become irrelevant, leading to the concept of 
universality classes for the critical exponents. Although first-order phase transitions are 
very common in nature, they have attracted much less attention from the theoretical 
point of view; in particular, the correlation length remains finite at such transitions, and 
so a universal behavior is not expected. Many interesting features occur in first-order 
phase transitions, like a discontinuity in the order parameter and the presence of latent 
heat; in addition to those, in some systems, it is possible for a first-order phase transition 
to become continuous due to changes in an external parameter, or to the presence of 

n 

microscopic inhomogeneities that may generate significant roundings in the transition . 
Therefore, one important question concerns the effects of random quenched impurities on 
the thermodynamic phase transition. 

The random-field Ising model (RFIM) was introduced by Imry and Ma [sj] and has 
attracted a lot of interest after the identification of its physical realizations. Probably 
the most important physical conception of the RFIM comes out to be a diluted Ising 
antiferromagnet in the presence of a uniform magnetic field 6|, LZl]. In these systems one 
has local variations in the sum of exchange couplings that connect a given site to other 
sites of the system, leading to local variations of the two-sublattice site magnetizations, 
and as a consequence, one may have local magnetizations that vary in both sign and 
magnitude. In these identifications, the effective random field at a given site is expressed 
either in terms of its local magnetization gI, or of the sum of the exchange couplings 
associated to this site 3]. 

Since then, many diluted antiferromagnets have been investigated, in such a way that 
systems like FexZni_xF2 and FexMgi_xCl2, for certain ranges of the concentration x (es- 



sentially high values of x, as mentionec 
experimental realizations of the RFIM 



^above) are nowadays considered as standard 
These systems are characterized by large 
crystal-field anisotropics, so that they may be reasonably well-described in terms of Ising 
variables. In particular, the compound FexMgi_xCl2 behaves like an Ising spin glass for 
X < 0.55, and is considered as a typical RFIM for higher magnetic concentrations. In the 
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RFIM regime, it presents a very curious behavior: one finds a first-order transition turn 



lOj; the concentration 



ing into a continuous one due to a change in the random fields 
at which the first-order transition disappears is estimated to be x = 0.6. 

From the theoretical point of view, many important questions remain open. At the 
mean-field level, it is well known that different probability distributions for the random 
fields may lead to distinct phase diagrams, e.g., a Gaussian probability distribution yields a 
continuous ferromagnetic-paramagnetic boundary {ifj], whereas for a bimodal distribution, 
this boundary exhibits a continuous piece at high temperatures ending up at a tricritical 



point, followed by a first-order phase transition at low temperatures 



12| . Indeed, it was 



argued that whenever an analytic symmetric distribution for the fields presents a minimum 



at zero field, one shou 
low temperatures 

r 

other authors 1 
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d expect a tricritical point and a first-order transition for sufficiently 
; such an argument has been considered further and improved by 



14l |. The presence of such a first-order transition for sufficiently low 



temperatures in short-range-interaction models, represents a point that has not been 
"uUy elucidated. For the three-dimensional RFIM, high-temperature series expansions 
15 1 and a zero-temperature scaling analysis IgI find continuous transitions for both 



Gaussian and bimodal distributions. Nevertheless, several zero-temperature studies of the 
Gaussian three-dimensional RFIM 17 , Q, [l^ suggest a first-order transition for such a 
model; this occurs because the magnetization and specific-heat critical exponents are very 
small, and so, it is difficult to determine whether the magnetization vanishes continuously, 
or discontinuously, at the transition. However, in four dimensions a zero-temperature 
analysis IgI leads to a first-order transition in the bimodal case and a continuous one for 
a Gaussian distribution, in agreement with the mean-field predictions. 

The crossover from first-order to continuous phase transitions has been investigated 



0,y, 



21 



22l |. One possible mechanism used 



through different theoretical approaches [id, 
to find such a crossover, or even to suppress the first-order transition completely, consists 
in introducing an additional kind of randomness in the system, e.g., bond randomness 



20 



2l| . By considering randomness in the field only, this crossover has 



Deen also ana- 



lyzed through zero-temperature studies, either wit 



simulations on a three-dimensional lattice 
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lin mean-field theory 22| , or numerical 



22 



. Since the effects observed in the di- 



luted antiferromagnet FexMgi_xCl2 occur for rather low temperatures (typically around 4 
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K) , they may be described satisfactorily by zero-temperature approaches. However, such 
zero-temperature approaches would not be appropriate for explaining similar effects at 
higher temperatures, that may possibly be found on similar systems. 

In the present work we present a RFIM that is able to exhibit a first-order phase 
transition turning into a continuous one, and characterized by a destruction of the first- 
order phase transition due a change in the random fields, at finite temperatures. The 
interactions among spins are of infinite range, i.e., in the limit where the mean-field 
approach is exact, and the random fields are defined through a double Gaussian probability 
distribution, which recovers in certain limits, the simple Gaussian and bimodal probability 
distributions. In the next section we present the model, find its free energy and phase 
diagrams. In section 3 we present our conclusions. 



II. THE MODEL AND ITS SOLUTION 



Let us define the infinite-range-interaction Ising model in the presence of an external 
random magnetic field, in terms of the Hamiltonian 

H ^ -—"^SiSj -^HiSi , (1) 



where Si = ±1 (i = 1,2, ...,iV) and the sum ^(jj) applies to all distinct pairs of spins. 
The random fields {Hi} are quenched variables, following a double Gaussian probability 
distribution, 

^1}. m 



pm 



exp 



[Hi 



HaY 



2a^ 



+ exp 



The probability distribution above depends on two parameters, Hq and cr, and modifies 
its form according to the ratio Hq/u, as exhibited in Fig. 1. Such a distribution is double- 
peaked for Hq/(t > 1, presents a single peak ior H^/a < 1, changing its concavity at the 
origin when Ho/a — 1. Besides that, in the limits Hq ^ and cr — > 0, one recovers the 
symmetric Gaussian and bimodal probability distributions, respectively. Changes between 
these two limits may be followed by analyzing the moments < if" >h {n = 1, 2, 3, • • • ) 
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FIG. 1: The probability distribution of Eq. ([2]) for typical values of the ratio Ho/a: (a) Ho/a = 
1/2,1,3/2; (b) Ho/a = 15. 

and in particular, through the kurtosis, k =< Hf >h /[3(< Hf >h)'^], that varies from 
K = 1/3 (bimodal limit) up to k = 1 (Gaussian limit), approaching unit in the limit 
Ho/cr — > 0, in which case one gets a perfect Gaussian distribution. For finite values of 
Ho /a one gets 1/3 < k < 1, and in particular, for the cases exhibited in Fig. 1 one 
has that n ^ 0.97 {Ho/ a = 1/2), n ^ 0.83 {Ho/cr = 1), k ^ 0.68 {Ho/(t = 3/2), and 
K ^ 0.34 {Ho/a = 15). 

For a given realization of the site fields {H^}, one has a corresponding free energy 
F{{Hi}), in such a way that the average over disorder, [F{{Hi})]H, becomes 



[F{m)]H = J ll[dH.P{Hd]F{{H,}) . 



(3) 



One can now make use of the replica method 23|, |2J] in order to get the free energy per 
spin as 



Pf = lim h\nZ{{H,})]H = lim Hm -^([^1// " 1) 

Af— >oo iV TV— >oo n— >0 iV n 



(4) 



where is the partition function of n copies of the original system defined in Eq. 
and (3 = l/{kT). Standart calculations lead to 



Pf = lim — min g{m") 



(5) 
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where 



^ - 1 In Tr„ exp(H+ ) - ^ In Tr„ exp(H^g) 

Q \ a / a 



(6) 



(7) 



In the equations above, the index a {a = 1, 2, n) is a rephca label and Tr^ represents a 
trace over the spin variables of each replica. The extrema of the functional g{m'^) yields 
the equilibrium equation for the magnetization of replica a, 



= i < >+ +i < >_ , 
2 +2 



where < >± refer to thermal averages with respect to the "effective Hamiltonians" Ti^^ 
in Eq. (El). 



If one assumes the replica-symmetry ansatz 23|,|2J], i.e., m" = m (V a), the free energy 
per spin [cf. Eqs. ©-([Tj)] and the equilibrium condition, Eq. ([8]), become 



/ 



J 2 11 

—m ^ , 

2 2p^J_ 

1 1 



dze ''^/^ In 2 cosh $" 



2/3v^y_ 



dze ^ In 2 cosh 



1 1 



m 



rfze"^'/^ tanh $+ + — ^ / dze~'^^^ tanh $" , 



(9) 
(10) 



where $^ = /5(Jm + a^; ± Hq). 

It is important to mention that t 
with the replica-symmetric solution 



le present system exhibits no instability associated 



25| . which usually appears due to p arameters char- 

23I, y. In the RFIM, 



acterized by two replica indices, like in the spin-glass problem 
one has a single phase transition associated with the magnetization and two phases are 
possible, namely, the ferromagnetic (m 7^ 0) and the paramagnetic (m = 0) ones. The 
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critical frontier separating these two phases may be found by solving Eq. ( fTOl) ; in the case 
of first-order phase transitions, we shall make use of the free energy per spin, Eq. (Q, as 
well. Let us then expand Eq. (JTOi) in powers of m, 



m 



Am + Bm^ + Cm^ + 0{m^) , (11) 



where the coefficients are given by 



A 
B 

C 



Ml -Pi} 



3 

15 



{l-4pi + 3p2} , 
{2 - 17pi + 30p2 - ISps} 



(12) 
(13) 

(14) 



with 



Pk 



1 r+co 

^= / (ize-"'/'tanh2'=/5(iJo + cT2) 



(15) 



The continuous critical frontier is determined by setting A = 1, provided that i? < 0. If 
a first-order critical frontier also occurs, the continuous line ends when 5 = 0; in such 
cases, the continuous and first-order critical frontiers meet at a tricritical point, whose 
coordinates may be obtained by solving the equations A = \ and 5 = numerically. In 
addition to that, for 5 > 0, the first-order critical frontier may be found by equating the 
free energies at each side of this line, i.e., /(m = 0) = /(m ^ 0). 

Using this procedure, we have calculated numerically the critical frontiers separating 
the paramagnetic and ferromagnetic phases for typical values of a/J (see Fig. 2). As will 
be seen below, the existence of a tricritical point at finite temperatures is restricted to 
the condition < (a/J) < ^J'l|e■n. At the threshold value (a/J) = ^J2|en = 0.4839 the 
tricritical point occurs at zero temperature. 

Let us now analyze how the above-mentioned critical frontier evolves along the zero- 
temperature axis; at T = 0, the free energy and magnetization become, respectively. 
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FIG. 2: Critical frontiers separating the paramagnetic (P) and ferromagnetic (F) phases of the 
infinite-range-interaction ferromagnet in the presence of a double Gaussian random field, for 
typical values of a/ J. Full (dashed) lines represent continuous (first-order) phase transitions; 
the black circles are tricritical points. Notice the colapse of the tricritical point at the zero- 
temperature axis for a/ J = 0.4839. 



/ 



J 2 

-— m 
2 



a 



Ho 
2 



exp 



erf 
( Jm 



Jm + Hn 



aV2 



2a2 



erf 



+ exp 



Jm — Ho 

(Jm - HoY 
2a2 



(16) 



1 JJm + Ho\ , 1 JJm-Ho 
m = -en ^ — H — eri ■= — 



(17) 



A procedure similar to the one described above for finite temperatures applies in this case, 
in such a way that one may expand Eq. (fT7|l in powers of m, 



m = am + bm^ + cm^ + 0{m'^) 



where 
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FIG. 3: Phase diagram of the infinite-range-interaction ferromagnet in the presence of a double 
Gaussian random field at zero temperature. The full (dashed) line represents a continuous 
(first-order) phase transition; the black circle denotes a tricritical point. 




A continuous critical frontier occurs at zero temperature for 6 < (i.e., Hq/ct < 1), in 
such a way that the condition a — 1 yields a relation involving Hq/J and a/ J for this 
critical frontier, 




One notices that the zero-temperature value of Hq/J decreases, for increasing values of 
(j/ J, and in particular, when a / J = \/2pn, one gets Hq/ J = as a solution. The critical 
frontier presents a tricritical point at zero-temperature given by 6 = 0, a = 1, 
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— = 1, - = W— = 0.4839. (23) 
a J \ en 



This represents the threshold value of a for the existence of a tricritical point; above this 
value one gets that Hq/ct < 1, in such a way that the probability distribution of Eq. ([2]) 
presents a single maximum at the origin (cf. Fig. 1) and there is no tricritical point, 
i.e., the critical frontier is completely continuous, in agreement with previous analyses 

For b > 0, one gets a first-order critical frontier at zero-temperature, which is associated 
with a tricritical point at finite temperatures. The phase diagram of the model, at zero 
temperature, is presented in Fig. 3. Similarly to what happened in the finite-temperature 
phase diagram, the first-order critical frontier changes into a continuous one, for increasing 



values of a/ J, like other zero-temperature studies of the RFIM jlOl. l22| 



III. CONCLUSIONS 

Summarizing, the main effects produced by an increase of a/ J in the present model are: 
(i) a decrease in the extension of the ferromagnetic phase; (ii) in the range of cr/ J for which 
there is a first-order transition line, one observes also a decrease in the extension of such a 
line; (iii) for sufficiently large values of a/ J, the first-order transitions are transformed into 
continuous ones. At the threshold value (cr/J) = y^2/e-K = 0.4839, the tricritical point 
where the continuous and the first-order transition lines meet, occurs at zero temperature, 
and for greater values of cr/J, the critical frontier is completely continuous, i.e., there is 
no tricritical point. Therefore, the ratio cr/J is directly related to the disorder in a real 
system; for the case of a diluted antiferromagnet, an increase in cr/J should play a similar 
role as an increase in the dilution. 

A crossover of the phase transition from first-order to a continuous one, due to an 
increase in the amount of disorder, has been observed in the diluted antiferromagnet 
FexMgi_xCl2, with 0.7 < x < 1.0 [ol, Q. Such an effect, that has been observed in low 
- but finite temperatures - has been explained in terms of zero-temperature analyses of 
different formulations of the RFIM [lol. I22I1. We believe that the present model is more 
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appropriate for a theoretical description of this effect. In this case, an increase in the 
measure of randomness in our model, cr/ J, would be related to a decrease in the magnetic 



concentration x, in such a way that the threshold value (o"/ J) = ^/ije-K would correspond 
to the critical value x = 0.6, at which the first-order transition disappears. In addition to 
that, the present model may also explain similar effects that could possibly be found, at 
higher temperatures, on other diluted antiferromagnets. 

Finaly, we argue that the double-Gaussian probability distribuition, defined above, is 
suitable for an appropriate theoretical description of the RFIM, being a better candidate 
for such a purpose than the two most commonly used distributions in the literature. 

(i) In the identifications of the RFIM with diluted antiferromagnets in the presence 
of a uniform magnetic field, the local random fields are expressed in terms of quantities 
that vary in both sign and magnitude 0, H]; this characteristic rules out the bimodal 
probability distribution from such a class of physical systems. 

(ii) Although the RFIM defined in terms of a simple Gaussian probability distribution 
for the fields is physically acceptable, it usually leads to a continuous phase transition 
at finite temperatures, either within mean- field 0, Q, 14 1, or standard short-range- 
interaction approaches 15|, [iGj. Such a system is not able to exhibit first-order phase 
transitions and tricritical points, that may occur in some diluted antiferromagnets [8|]. 

(iii) By varying appropriately the ratio Ho/a (a ratio related to the external applied 
uniform field and the dilution in a real system) in the double-Gaussian probability dis- 
tribution of the present RFIM, one may adjust the model to given physical situations, in 
order to reproduce a wide variety of physical effects that occur in diluted antiferromagnets, 
like continuous and first-order phase transitions, as well as tricritical points. 
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